ARL-TR-94-21 


Copy  No. 


Bottom  Backscatter  from  Trapped  Bubbles  -  II 

Final  Report  under  Contract  N00039-91-C-0082 
TD  No.  01A2063,  Bottom  Backscatter  from  Treipped  Bubbles  -  II 


Frank  A.  Boyle 
Nicholas  P.  Chotiros 


Applied  Research  Laboratories 

The  University  of  Texas  at  Austin 
P.  O.  Box  8029  Austin,  TX  7871 3-8029 


1  December  1 994 
Final  Report 

15  April  1993-15  April  1994 

Approved  for  public  release; 
distribution  is  unlimited. 


Prepared  for: 

Naval  Research  Laboratory 
Stennis  Space  Center,  MS  39529-5004 


Monitored  by: 

Space  and  Naval  Warfare  Systems  Command 
Department  of  the  Navy 
Arlington,  VA  22245-5200 


19950712  001 


pn® 


S 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


the^llection  of  Information.  Send  comments  regarding  this  burden  estimate 
is  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
>-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction 


Project  (0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 


1  Dec  94 


4.  TITLE  AND  SUBTITLE 


Bottom  Backscatter  f rom  Trapped  Bubbles  -  II,  Final  Report  under  Contract 
N00039-91-C-0082,  TD  No.  01A2063,  Bottom  Backscatter  from  Trapped  Bubbles  ■  II 


S.  FUNDINO  NUMBERS 


Contract  N00039-91  -C-0082 
TD  No.  01A2063 


6.  AUTHOR(S) 

Boyle,  Frank  A. 


Chotiros,  Nicholas  P. 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  AODRESS(ES) 

Applied  Research  Laboratories 
The  University  of  Texas  at  Austin 
P.O.  Box  8029 
Austin,  TX  78713-8029 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-94-21 


9.  SPONSORING/MONITORINQ  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory  Space  and  Naval  Warfare 


Naval  Research  Laboratory 
Stennis  Space  Center,  MS  39529-5004 


Systems  Command 
Department  of  the  Navy 
Arlington,  VA  22245-5200 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

A  model  for  acoustic  backscatter  from  trapped  gas  bubbles  in  sandy  sediments  was  described  in  "Bottom  Backscatter 
from  Trapped  Bubbles,"  Applied  Research  Laboratories  Technical  Report  No.  93-15  (ARL-TR-93-15).  In  that  model, 
trapped  bubbles  were  assumed  to  scatter  as  if  they  were  free  bubbles  in  open  water.  In  this  report,  the  effects  of  bubble 
confinement  in  sediment  pores  on  the  resonance  behavior  of  the  bubble  are  accounted  for.  This  is  done  by  assigning 
the  pore  fluid  an  effective  density  that  differs  from  its  actual  density,  accounting  for  the  fact  that  the  fluid  is  partially 
confined  within  pores.  The  effective  density  is  computed  by  way  of  the  Biot  theory.  Two  effective  densities  are 
specified,  one  for  each  of  the  two  compressional  waves  that  the  Biot  theory  predicts.  As  a  result,  the  medium  has  two 
scattering  cross  sections,  which  are  both  included  in  the  resulting  expression  for  scattering  strength. 


14.  SUBJECT  TERMS 

acoustic  impedance 
backscatter  ( 

Biot 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


15.  NUMBER  OF  PAGES 


bubble 

effective  density 


pore  pressure 
reciprocity 
sandy  sediments 


scattering  cross  section  _ 29_ 

shallow  grazing  angle  16.  price  code 


18.  SECURITY  classification 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Sid.  239-18 
298-102 


This  page  intentionally  left  blank. 


TABLE  OF  CONTENTS 


Page 


LIST  OF  FIGURES .  v 

PREFACE .  vii 

EXECUTIVE  SUMMARY .  ix 

1.  INTRODUCTION  . 1 

2.  TRAPPED  BUBBLE  BACKSCATTERING  STRENGTH  .  3 

2. 1  BACKSCATTERING  CROSS  SECTION  OF  A 

CONSTRAINED  BUBBLE .  3 

2.2  RELATIONSHIP  BETWEEN  INCIDENT  AND 

BACKSCATTERED  PRESSURE .  5 

2.2.1  Relationship  between  Acoustic 

Pressures  and  Fluid  Velocities .  5 

2.2.2  Acoustic  Impedance .  7 

2.2.3  Backscattered  Pressure  from  an  Element 

dxdydz  of  Sediment  Volume  by  Reciprocity .  7 

2.3  CALCULATION  OF  PORE  PRESSURES  AND 

EFFECTIVE  DENSITIES  WITH  THE  BIOT  THEORY .  1 0 

2.4  SEDIMENT  INTERFACE  BACKSCATTERING  STRENGTH .  12 

3.  COMPARISON  OF  THEORY  WITH  EXPERIMENT .  17 

4.  SUMMARY  AND  CONCLUSIONS  .  21 

REFERENCES . 23 


Accesion  For 

NTIS  CRA&I 
OTIC  TAB 
Unannounced 
Justification 


By . . . 

Distribution  / 


Availability  Codes 

Dist 

fi-l 

Avail  ; 
Spt 

3rid/or 

cial 

This  page  Intentionally  left  blank. 


IV 


LIST  OF  FIGURES 


Figure  Pas© 

2.1  Calculation  of  backscattered  pressure  by  reciprocity .  8 

2.2  Backscatter  from  the  volume  below  surface  element  dxdy .  1 3 

2.3  Incident  and  scattered  pressures  as  defined  for  an 

interface  scattering  element . .  1 5 

3.1  Comparison  of  experimental  data  with  backscattering 

strengths  predicted  by  Eq.  (2.36) .  1 8 


V 


This  page  intentionally  left  blank. 


VI 


PREFACE 


This  document  is  the  final  report  on  work  that  Applied  Research 
Laboratories,  The  University  of  Texas  at  Austin  (ARLUT),  was  tasked  to  perform 
under  Contract  N00039-91-C-0082,  TD  No.  01A2063,  Bottom  Backscatter  from 
Trapped  Bubbles  -  II. 


vii 


This  page  intentionally  left  blank. 


EXECUTIVE  SUMMARY 


In  "Bottom  Backscatter  from  Trapped  Bubbles"  by  Boyle  and  Chotiros^  a 
model  for  shallow  grazing  angle  acoustic  backscatter  from  sandy  sediments 
was  developed.  It  included  a  Biot  model  for  acoustic  propagation  and  a 
scattering  mechanism  from  trapped  bubbles  in  the  sediment  pores.  One  of  the 
assumptions  made  was  that  the  resonance  frequencies  and  scattering  cross 
sections  of  trapped  bubbles  are  the  same  as  those  the  bubbles  would  have  if 
they  were  surrounded  by  open  water.  The  possible  effects  of  confinement 
within  sediment  pores  were  neglected.  In  this  report,  a  revised  expression  for 
the  sediment  backscattering  strength  is  developed  that  includes  the  effects  of 
bubble  confinement  in  pores. 

The  Biot  slow  and  fast  compressional  waves  are  considered  separately. 
Shear  waves  are  neglected  because  they  do  not  couple  strongly  into  the  pore 
fluid.  There  are  separate  acoustic  impedances  and  scattering  cross  sections, 
specific  for  the  fast  and  slow  waves.  The  scattering  problem  is  treated 

separately  for  each  of  the  two  wave  types.  The  total  scattering  cross  section 
includ©s  compon©nts  from  both  wav©s. 

The  model  requires  a  number  of  input  parameters.  There  are  methods  to 
measure  or  estimate  all  but  two.  They  are  (1)  the  sediment  pore  fluid  gas 

fraction  and  (2)  the  bubble  radius  to  pore  radius  ratio  (rbp).  These  are  treated  as 
free  parameters. 

Model  predictions  are  compared  with  backscatter  data  from  four  sites. 
The  model,  by  adjustment  of  its  free  parameters,  can  be  made  to  fit  the 
observations.  This  supports  the  hypothesis  that  trapped  gas  bubbles  may  be  a 
significant  source  of  acoustic  backscatter  in  sediments. 

Very  small  gas  fractions  of  10-®  to  10-6  are  sufficient  to  fit  observed  levels 
of  backscatter.  There  is  no  practical  way  to  measure  such  small  gas  fractions  at 
present,  so  an  independent  verification  is  not  presently  available. 


The  bubble  size  to  pore  size  ratio  required  to  fit  grain  size  dependence  of 
existing  backscatter  data  is  unreasonably  large.  This  ratio  also  increases 


IX 


steadily  as  the  sediment  grain  size  decreases.  This  behavior  casts  doubt  on  the 
accuracy  of  the  current  method  of  estimating  trapped  bubble  sizes,  specifically, 
that  the  bubble  size  distribution  is  proportional  to  the  pore  size  distribution.  A 
more  accurate  way  of  estimating  bubble  size  distribution  is  needed. 

The  values  of  the  input  parameters  that  fit  the  model  to  existing 
backscatter  measurements  appear  to  be  reasonable.  The  fact  that  the  bubble  to 
pore  radius  ratio  is  greater  than  unity  suggests  that  other  factors  besides  pore 
size  influence  the  size  of  a  trapped  bubble.  When  a  more  accurate  estimate  of 
the  bubble  size  distribution  becomes  available,  it  can  be  incorporated  into  the 
model.  Until  then,  the  current  technique  allows  a  semi-empirical  estimate  of 
bubble  size  distribution,  which  results  in  a  reasonable  fit  of  the  model  to 
observed  backscatter  data. 


X 


1.  INTRODUCTION 


The  objective  of  this  work  is  an  improved  model  for  high  frequency 
acoustic  backscatter  from  sandy  sediments  at  shallow  grazing  angles.  Initial 
development  of  the  model  was  described  by  Boyle  and  Chotiros  in  1 992.^  In  this 
initial  work,  a  comprehensive  compilation  of  shallow  grazing  angle  backscatter 
measurements  was  presented.  Some  of  the  data  featured  a  broad  maximum  in 
the  backscattering  strength  spectrum  that  was  not  explainable  with  current 
models.  A  new  hypothetical  mechanism  for  acoustic  backscatter  was  proposed, 
involving  resonance  scattering  from  trapped  gas  bubbles.  A  preliminary 
backscatter  model  based  on  this  scattering  mechanism  was  then  developed.  In 
1993,  a  more  complete  mode|2  was  developed  that  combined  a  Biot  propagation 
model  with  a  trapped  bubble  backscatter  model.  The  model  included  a  means  of 
estimating  the  trapped  bubble  size  distribution  from  the  grain  size  distribution, 
based  on  the  assumption  that  the  bubble  size  distribution  mirrors  the  pore  size 
distribution. 

One  of  the  assumptions  made  in  the  1993  model  is  that  trapped  bubbles 
respond  acoustically  as  if  they  are  free  bubbles  surrounded  by  water.  In  this 
report,  the  acoustic  response  of  bubbles  that  are  constrained  within  sediment 
pores  is  investigated.  The  presence  of  solid  particles  in  the  fluid  surrounding  the 
bubbles  is  accounted  for.  The  result  is  a  new  model  with  no  assumptions 
regarding  the  bubbles'  confinement  in  the  sediment. 

This  report  is  arranged  as  follows:  in  Sec.  2  a  theoretical  model  is 
developed  to  describe  the  acoustic  response  of  bubbles  trapped  in  sediment. 
Section  3  contains  comparisons  between  model  predictions  and  experimental 
measurements.  A  discussion  of  the  results  follows  in  Sec.  4. 
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2.  TRAPPED  BUBBLE  BACKSCATTERING  STRENGTH 


The  objective  of  the  backscatter  model  is  an  expression  for  the 
backscattering  strength  of  the  sediment  interface,  due  to  scattering  from  bubbles 
in  the  volume  below  the  interface.  The  acoustic  response  of  the  bubbles  is 
influenced  by  the  fact  that  the  bubbles  are  confined  in  pores  between  sand 
grains. 


Sediment  particles  influence  the  acoustic  behavior  of  bubbles  because 
they  affect  the  effective  density  of  the  surrounding  fluid  medium.  This  effective 
density  is  the  density  of  an  equivalent  unbounded  fluid,  in  which  bubbles  would 
have  the  same  acoustic  response  that  they  have  in  sandy  sediment. 

The  influence  of  the  effective  density  manifests  itself  in  two  ways.  First, 
the  backscattering  cross  section  of  each  bubble  is  affected  by  the  density  of  the 
surrounding  fluid.  Second,  the  relationship  between  incident  and  backscattered 
pressure,  which  is  developed  from  the  reciprocity  principle, 3-^  is  affected  by  the 
impedance  difference  across  the  sediment  boundary.  The  two  effects  are 
discussed  separately  in  Secs.  2.1  and  2.2  below.  Section  2.3  is  a  discussion  of 
how  pore  fluid  effective  densities  and  pressures  are  determined  by  way  of  the 
Biot  theory.  In  Sec.  2.4  the  pore  pressure  is  used  to  obtain  an  expression  for  the 
sediment  backscattering  strength. 

In  the  following  analysis,  complex  quantities  are  printed  in  boldface.  Real 
quantities  are  plain  text. 

2.1  BACKSCATTERING  CROSS  SECTION  OF  A  CONSTRAINED  BUBBLE 

In  this  section,  an  expression  is  presented  for  the  acoustic  scattering  cross 
section  of  a  bubble  in  terms  of  the  density  of  the  surrounding  fluid.  An 
expression  for  the  appropriate  value  for  this  density  is  later  presented  in  Sec.  2.3. 

The  backscatter  model  described  in  Ref.  2  uses  the  expression  for  the 
scattering  cross  section  of  a  single  spherical  bubble,  given  by  Wildt^: 


3 


Or(rb)  = 


(2.1) 


where  Cr  is  the  scattering  cross  section,  rb  is  the  bubble  radius,  f  is  the  incident 
acoustic  frequency,  fr  is  the  resonance  frequency,  and  8  is  the  damping  constant, 
given  by 


5  =  krb  +  ^T 


P^r^b 


where  k  is  the  acoustic  wavenumber  in  the  surrounding  medium;  p  is  the  shear 
viscosity  of  the  fluid;  p  is  the  density  of  the  medium  surrounding  the  bubbles; 
ci)r=27cfr  is  the  angular  resonance  frequency;  d/b  is  given  by 


d 

b 


=  3(r-i) 


x{sinhX  +  sinX]-2(coshX  -  cosxj 
X^(coshX  -  cosxj  +  3(y-1  )x|sinhX  -  sinX] 


(2.3) 


where  y  is  the  ratio  Cp/Cv  of  specific  heats  for  the  gas  inside  the  bubble;  and  X  is 
given  by 


X  =  rb 


(2.4) 


The  density  of  the  gas  inside  the  bubble  is  pg,  Cpg  is  the  specific  heat  of  the  gas 
at  constant  pressure,  and  Kg  is  the  thermal  conductivity  of  the  gas. 


The  resonance  frequency  fr  is  given  by® 


(2.5) 
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where  Pq  is  the  ambient  static  pressure  and  p  is  the  density  of  the  surrounding 
medium:  b  and  p  are  quantities  which  must  be  included  to  account  for  the 
surface  tension  x  of  the  bubble  wall  and  the  thermal  conductivity  Kg  of  the  gas 
inside  the  bubble: 


b  = 


X 


f  sinhX-sinX 
VcoshX-cosX 


(2.6) 


p  =  1+ 


(2.7) 


The  ambient  density  p  appears  in  Eqs.  (2.2)  and  (2.5).  In  the  simple  case 
of  a  pure  fluid  medium,  p  is  simply  the  density  of  the  fluid.  When  the  medium 
surrounding  the  bubble  contains  sediment  particles,  the  effective  fluid  density 
presented  in  Sec.  2.3  must  be  used. 

2.2  RELATIONSHIP  BETWEEN  INCIDENT  AND  BACKSCATTERED 
PRESSURE 


In  this  section,  a  general  expression  for  the  backscattered  pressure  from 
an  element  of  sediment  volume  is  derived.  In  Sec.  2.2.1  the  relationship  between 
pore  fluid  acoustic  pressure  and  velocity  is  presented.  In  Sec.  2.2.2,  an 
expression  for  the  pore  fluid  acoustic  impedance  is  developed.  In  Sec.  2.2.3  the 
principle  of  acoustic  reciprocity  is  used  to  develop  an  expression  for  the 
backscattered  pressure  from  a  sediment  volume  element.  This  backscattered 
pressure  is  expressed  in  terms  of  the  pressure  incident  upon  the  scattering 
element,  which  can  be  obtained  with  the  Biot  model,  described  in  Sec.  2.3. 

2.2.1  Relationship  between  Acoustic  Pressures  and  Fluid  Velocities 


The  linear  wave  equation  in  a  fluid  medium  is  given  by 


v2p  = 


(2.8) 
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where  c  is  the  phase  velocity  of  a  compressionai  wave  in  the  fluid.  For  a 
pressure  field  that  is  sphericaiiy  symmetric  about  a  source  at  the  origin,  there  is 

only  radial  dependence.  The  operator  in  this  case  is  given  by 

=  ■  (2-9) 

Upon  substitution  of  Eq.  (2.9)  into  Eq.  (2.8)  and  rearrangement,  the  wave 
equation  can  be  written  in  the  form 


a^(rp)  _  1  3^(rp) 

ar2  "c2  at2 


(2.10) 


Generai  solutions  for  the  quantity  rp  can  be  expressed  as 


rp  =  fi(ct-r)  +  f2(ct+r) 


(2.11) 


where  the  first  term  represents  a  general  incoming  wave  and  the  second  term  a 
general  outgoing  wave.  If  we  consider  outgoing  harmonic  waves  only,  the 
pressure  field  can  be  expressed  as 


p  =  yexp(i(kr-cot)) 


(2.12) 


where  A  is  an  arbitrary  complex  amplitude,  k=(o/c+ia  is  the  acoustic 
wavenumber,  a  is  the  absorption  in  Np/m,  and  (o=2ni  is  the  angular  frequency. 
For  small  amplitude  signals  the  fluid  velocity  is  related  to  the  acoustic  pressure 
according  to  Euler's  equation, 7 

Po^  =  “^P  >  (2-13) 

where  po  is  the  average  density  of  the  medium.  Upon  combination  of  Eqs.  (2.12) 
and  (2.13),  a  relation  between  acoustic  pressure  and  fluid  velocity  is  obtained. 
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kr)^ 


(2.14) 


where  f  is  the  unit  vector  in  the  radial  direction. 


2.2.2  Acoustic  Impedance 


The  acoustic  impedance  is  defined  in  terms  of  the  acoustic  pressure  and 
fluid  velocity, 

z,  =  ^  ,  (2.15) 


where  V|  is  the  velocity  component  in  the  direction  of  the  velocity  vector  v  .  For 
small  amplitude  signals  Eq.  (2.15)  simplifies  to 


Zi  =  — 


V| 


(2.16) 


Equations  (2.14)  and  (2.16)  can  be  combined  to  form  an  expression  for  the 
acoustic  impedance: 


7  _  Pf 

4-  _  - - - 


K) 


(2.17) 


At  large  distances  from  the  source,  where  kjr  is  large,  the  impedance  of 
Eq.  (2.17)  reduces  to  the  plane  wave  acoustic  impedance  Z\  =  pjC. 

2.2.3  Backscattered  Pressure  from  an  Element  dxdydz  of  Sediment 
Volume  by  Reciprocity 

Consider  the  situation  where  the  scatterer  is  an  element  of  sediment 
volume  dxdydz.  In  this  case  it  is  convenient  to  surround  the  source  and  scatterer 
with  virtual  spheres  with  radius  r  much  greater  than  the  acoustic  wavelength  X,  as 
illustrated  in  Fig.  2.1. 
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Source 


Water 


(a)  Acoustic  propagation  of  incident  sound. 


Water 


(b)  Acoustic  propagation  of  backscattered  sound. 


Figure  2.1 

Calculation  of  backscattered  pressure  by  reciprocity. 

Reciprocity  is  applied  to  virtual  spheres  of  equal  radius  about  source  and  scatterer.  The  radius 
is  much  greater  than  the  acoustic  wavelength  so  that  plane  wave  acoustic  impedances  can  be 
used,  and  much  smaller  than  the  separation  between  source  and  scatterer  so  that  the  projected 
pressures  Ip^l  and  IP2I  are  approximately  constant  across  the  spheres'  surfaces. 
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The  spherical  surface  surrounding  the  source  has  a  surface  velocity  vo, 
which  generates  a  pressure  pi  at  the  location  of  the  scattered  The  scatterer 
responds  with  a  scattered  surface  velocity  vi,  which  induces  a  backscattered 
pressure  P2  at  the  source.  These  pressures  and  surface  velocities  can  be  related 
by  acoustic  reciprocity,  which  states  that,  in  a  linear  medium,  a  source  and 
receiver  can  be  swapped  with  no  change  in  the  ratio  of  received  to  transmitted 
signals.  This  swapping  of  positions  can  be  interpreted  to  represent  the 
backscatter  case,  where  the  scatterer  acts  as  a  projector  and  the  source  as  a 
receiver.  In  terms  of  pressures  and  velocities,  the  reciprocity  principle  states: 


PI-P2 

Vo" 


(2.18) 


By  inversion  of  Eq.  (2.16),  the  fluid  velocities  at  the  virtual  surfaces  of  the 
projector  and  scatterer  can  be  expressed  in  terms  of  local  pressures  and 
impedances: 


and 


Vo  = 


Vi  = 


Po 

Zo 

(2.19) 

ePi 

Zi  ’ 

(2.20) 

where  vo,  Po,  and  Zq  are  the  fluid  velocity,  the  acoustic  pressure,  and  the 
impedance  at  the  surface  of  the  virtual  sphere  surrounding  the  projector.  vi  and 
Zi  are  the  scattered  fluid  velocity  and  the  acoustic  impedance  at  the  virtual 
sphere  surrounding  the  scatterer.  e  is  a  transfer  function,  from  incident  pressure 
Pi  to  scattered  pressure  epi  at  the  surface  of  the  scatterer's  virtual  sphere.  It 
includes  the  effects  of  sound  pressure  generation  by  the  scattering  element 
dxdydz  as  an  effective  source,  as  well  as  propagation  of  this  pressure  to  the 
virtual  sphere  surface.  As  illustrated  in  Fig.  2.1 ,  this  propagation  must  take  place 
as  if  it  were  in  the  water  column.  If  we  neglect  the  attenuation  in  the  water 
column,  the  average  square  magnitude  of  e  is 


/|^|2\_<ybv  dxdydz 
4itri 


(2.21) 
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where  dxdydz  is  the  sediment  volume  element  and  rs  is  the  radius  of  the  virtual 
spheres  surrounding  source  and  scatterer.  Obv  is  the  volume  scattering  cross 
section,  defined  as  ratio  of  scattered  power  to  incident  intensity  for  the  sediment 
volume  scattering  element,  averaged  over  the  complete  ensemble  of  possible 
bubble  distributions  within  dxdydz.  By  combination  of  Eqs.  (2.18),  (2.19),  (2.20), 
and  (2.21),  an  expression  for  the  square  magnitude  of  the  acoustic  pressure 
returned  to  the  projector  from  the  scattering  element  dxdydz  is  obtained; 


<IP2R  =  | 


f 

Pi 

2  0bv 

Po 

47trf 

dxdydz 


(2.22) 


This  quantity  is  proportional  to  the  power  returned  to  the  projector  from  the 
scattering  element. 

Since  the  radius  of  the  virtual  spheres  was  assumed  large  in  comparison  to  the 
wavelength  X,  plane  wave  acoustic  impedances  can  be  used  for  Zo  and  Zi: 


|Zol 


Po 

Vo 


=  Po^o 


(2.23) 


and 


=  PlCi 


(2.24) 


where  po,  vo,  pi,  and  vi  are  plane  wave  acoustic  pressures  and  fluid  velocities  in 
the  water  column  and  sediment,  respectively,  po  and  pi  are  densities  of  the 
water  and  sediment  and  Co  and  ci  are  corresponding  phase  velocities. 

2.3  CALCULATION  OF  PORE  PRESSURES  AND  EFFECTIVE  DENSITIES 
WITH  THE  BIOT  THEORY 

For  a  sandy  sediment,  pi  and  vi  can  be  obtained  via  the  Biot  model, 
described  in  Ref.  2.  The  Biot  model  predicts  fast  and  slow  compressional  waves 
in  the  pore  fluid.  Since  shear  waves  do  not  couple  strongly  into  the  pore  fluid, 
they  can  be  neglected.  Each  wave  will  have  its  own  fluid  velocity  and  pressure. 
The  pore  fluid  is  therefore  assigned  two  partial  acoustic  impedances,^  one  for 
each  of  the  Biot  compressional  waves. 
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Zi,  =  ^  (2.25) 

and 

Zl.  =  .  (2.26) 

where  the  subscript  1  corresponds  to  the  pore  fluid  and  the  subscripts  f  and  s 
correspond  to  the  fast  and  slow  waves,  respectively. 

The  pore  fluid  also  has  two  effective  densities,  one  for  each  of  the  Biot 
waves  that  the  medium  supports; 

Pfast  “  (2-27) 

and 

Psiow~^~o^  ’  (2.28) 

where  Cf  and  Cs  are  the  fast  and  slow  wavespeeds,  obtainable  via  the  Biot  model. 

The  simultaneous  presence  of  the  fast  and  slow  acoustic  waves  affects 
the  form  of  the  returned  pressure  from  a  scattering  element,  Eq.  (2.22).  For 
scattering  within  a  Biot  medium,  the  expression  for  the  magnitude  squared 
backscattered  pressure  is 


dxdydz 


(2.29) 


where  pi  in  this  expression  is  the  total  pore  pressure,  Pif+pis,  due  to  fast  and 
slow  waves. 


^bvf  °^bvs  ensemble  average  volume  scattering  cross  se(^ion 

densities  for  fast  and  slow  waves,  respectively.  They  represent  the  collective 
backscatter  contributions  from  all  the  trapped  bubbles  in  the  sediment  and  have 
sensitive  dependence  on  the  size  distribution  of  trapped  bubbles.  The 
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distribution  of  bubble  radii  is  assumed  proportional  to  the  distribution  of  pore  radii, 
with  the  proportionality  factor  rbp: 

rbp  =  bubble  radius/pore  radius  .  (2.30) 

The  distribution  of  pore  radii  is  related  to  the  distribution  of  grain  radii,  which  can 
be  measured.  This  relationship  is  detailed  in  Ref.  2. 

An  expression  for  the  bubble  size  distribution,  on  which  these  scattering 
cross  sections  are  based,  is  detailed  in  Ref.  2  and  is  based  on  a  bubble  size 
distribution  that  mirrors  the  sediment  pore  size  distribution.  Given  a  pore  size 
distribution,  which  can  be  estimated  from  the  grain  size  distribution,  the  bubble 
size  distribution  is  specified  by  a  single  parameter,  rbp,  that  specifies  the  bubble 
size  to  pore  size  ratio. 

2.4  SEDIMENT  INTERFACE  BACKSCATTERING  STRENGTH 

In  this  section,  the  backscatter  contribution  from  a  sediment  volume 
element,  given  by  Eq.  (2.29),  is  used  to  determine  the  backscattering  strength  of 
the  sediment  interface.  In  order  to  specify  an  interface  backscattering  strength, 
scattering  from  within  the  sediment  volume  must  be  recast  as  "effective"  interface 
scattering.  As  illustrated  in  Fig.  2.2,  one  can  define  the  backscattered  intensity 
from  an  interface  element,  dxdy,  as  the  sum  of  backscattered  intensities  from  all 
volume  elements  dxdydz  below  dxdy. 

This  definition  of  the  interface  backscattering  contribution  is  based  on  the 
assumption  of  plane  waves  incident  upon  the  interface.  The  resulting 
backscatter  model  will  therefore  be  most  applicable  when  the  sediment  interface 
is  far  from  the  projector,  i.e.,  when  wavefront  curvature  is  negligible  in 
comparison  with  acoustic  penetration  depth  into  the  sediment. 

The  average  magnitude  squared  of  the  backscattered  pressure  Ipbsl^ 
returned  to  the  projector  from  the  sediment,  per  unit  area  of  interface,  is  given  by 
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Projector; 

backscattered  p  =  pbs 


Surface  element  dxdy 


Figure  2.2 

Backscatter  from  the  volume  below  surface  element  dxdy. 
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(|Pbsl^)=  pj 


Pi  1^  ^bvf  ^Oi 


Po|  1  47crf  |Zif' 


.  *^bvs  .^0. 

*  4nrl 


(2.31) 


The  backscattering  strength  of  an  interface  element  is  defined  as 

BS  =  10log^^^  , 

I  Pine  I 


(2.32) 


where  ps  is  the  backscattered  pressure  at  unit  distance  rim  from  the  interface 
element  and  pinc  is  the  pressure  incident  upon  interface,  as  illustrated  in  Fig.  2.3. 
These  pressures  are  related  to  Pbs  3rid  po  in  the  following  way: 


P,|=|Pbs|,'^e'“ 


(2.33) 


|Pinc|  =  |Po|7e 


s  ^-ra 


(2.34) 


where  r  is  the  separation  between  the  projector  and  the  interface  element  and  a 
is  the  absorption  coefficient  of  the  water  column  in  units  of  nepers  per  unit 
distance.  By  combining  Eqs.  (2.31),  (2.32),  (2.33),  and  (2.34),  the  backscattering 
strength  of  an  interface  element  of  unit  area  is  expressible  in  terms  of  po  and  pi: 


BS  =  10  log 


1*2  |.4 
*1m*s 


2  2 1 
|2  Obvf  Zq  *^bvs  .?0_n 

I  -W  *  4n  , 


+  40  log(e)  r  a 


(2.35) 


In  terms  of  the  incident  pressure  pinc  at  the  interface,  the  backscattering 
strength  is 
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Figure  2.3 

Incident  and  scattered  pressures 
as  defined  for  an  interface  scattering  element. 
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3.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 


Figure  3.1  is  a  comparison  of  model  predictions  with  backscatter 
measurements  from  four  selected  sites.9'^2  Backscattering  strengths  are  plotted 
against  normalized  grain  size,  which  is  proportional  to  frequency  as  follows: 


normalized  grain  size  = 


frequency  *  mean  grain  diameter 
sound  speed  in  water 


(3-1) 


The  input  parameters  for  each  site  are  listed  in  Table  3.1 .  In  each  case  a  best  fit 
was  obtained  by  varying  the  gas  fraction  ^  and  the  bubble  to  pore  radius  ratio 
(rbp).  The  resulting  gas  fractions  varied  between  1x10-®  and  1x10-®,  and  rbpS 
varied  between  1 .71  and  8.43. 


The  fact  that  the  bubble  to  pore  radius  ratios  rbp  were  greater  than  unity 
appears  to  have  the  unphysical  implication  that  trapped  bubbles  are  larger  than 
their  surrounding  pores.  One  possible  explanation  for  this  is  that  the  assumption 
that  bubble  size  distributions  mirror  the  pore  size  distributions  may  be  inaccurate. 
Surface  tension,  for  example,  may  force  very  small  bubbles  into  solution,  skewing 
the  distribution  toward  bubbles  of  larger  radii.  Such  an  effect  would  affect  bubble 
distributions  in  fine  grained  sands  more  severely  than  those  in  coarse  grained 
sands.  The  idea  that  the  bubble  size  distribution  may  favor  larger  radii  than  the 
pore  size  distribution  is  supported  by  the  fact  that,  as  shown  in  Table  3.1,  rbp 
increases  steadily  with  decreasing  mean  grain  size.  Other  possible  explanations 
for  the  large  rbpS  are  that  bubbles  in  neighboring  pores  coalesce  or  couple 
acoustically,  or  that  other  mechanisms  besides  trapped  bubbles  contribute  to  the 
backscattering  strengths. 

By  varying  the  two  free  parameters,  the  behavior  of  the  model  can  be 
adjusted  to  match  observations  well.  The  observed  peak  in  backscattering 
strength  is  matched  by  the  model,  given  appropriate  values  for  ^  and  rbp.  As  a 
result,  very  small  gas  fractions  are  inferred.  At  present  the  authors  are  aware  of 
no  practical  way  of  measuring  such  small  gas  fractions  in  sediments. 
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Backscattering  strength  -  dB  Backscattering  strength  -  dB 


Kings  Bay,  1985 

Gas  fraction . 1x10'® 

Bubble:pore  radius  ratio  (r^p) . 3.11 


Panama  City,  1986 
Gas  Fraction  IxlO"^ 


rbp . 8.43 


Jacksonville,  1989 
Gas  Fraction  2x10*® 
rbp . 1-71 


McKinney  and 
Anderson,  1964 
Gas  Fraction  1x1  O'® 
rbp .  8.43 


Figure  3.1 

Comparison  of  experimental  data  with 
backscattering  strengths  predicted  by  Eq.  (2.36). 
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4.  SUMMARY  AND  CONCLUSIONS 


A  model  has  been  developed  to  predict  acoustic  backscatter  from  gas 
bubbles  trapped  in  the  pores  of  a  sandy  sediment.  The  model  incorporates  a 
Biot  acoustic  propagation  model  and  allows  for  scattering  of  both  fast  and  slow 
waves.  The  effects  of  bubble  confinement  within  sediment  pores  on  bubble 
resonances  and  acoustic  scattering  cross  sections  are  included. 

The  new  model  has  two  free  input  parameters,  the  gas  fraction  and  the 
bubble  to  pore  radius  ratio.  Since  there  is  currently  no  method  of  measuring 
either  of  them,  they  are  simply  adjusted  to  fit  the  model  to  measured  data.  By 
adjusting  these  parameters,  it  is  possible  to  get  the  model  to  fit  data  collected 
experimentally  at  several  sites.  The  fit  includes  a  broad  peak  that  sometimes 
occurs  in  the  backscattering  strength  spectrum  that  cannot  be  explained  with 
other  models. 

A  very  small  gas  fraction  of  between  1 0'^  and  1 0  ®  was  sufficient  to  get  the 
model  to  match  the  observed  backscatter.  These  gas  fractions  are  smaller  than 
can  be  measured  at  present  so  there  is  no  independent  verification  of  the  model. 

There  are  reasons  to  expect  gas  in  shallow  water  sediments.  13  Gas 
bubbles  have  been  observed  in  core  samples  of  marine  sediments,  as  well  as  in 
echo  soundings,  which  have  recorded  bubbles  escaping  to  the  surface.  This 
suggests  that  trapped  bubbles  may  be  significant  in  high  frequency  acoustic 
backscatter  from  sediments. 

The  large  values  of  the  bubble  to  pore  radius  ratio  suggest  that  the 
assumed  shape  of  the  bubble  size  distribution  function  is  inaccurate.  The  fact 
that  rbp  increases  steadily  with  decreasing  mean  grain  size  suggests  that  the 
actual  bubble  size  distribution  contains  fewer  small  bubbles  and  more  large 
bubbles.  A  study  of  other  factors,  besides  pore  size,  that  might  affect  bubble  size 
may  allow  a  better  estimate  of  the  bubble  size  distribution. 
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